Nanoparticle networks as chemoselective sensing devices 
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We theoretically analyzed transport properties of a molecular network constructed of gold 
nanoparticles linked with oligophenylenevinulene (OPV) molecules. We showed that the conduc- 
tance of such system was strongly reduced when trinitrotoluene (TNT) became attached to the 
OPV linkers in the network. The reported results are based on the ab initio electronic structure 
calculations. These results corroborate and elucidate experiments which revealed significant drops 
in the conductance the network while the latter was exposed to TNT vapors. The results suggest 
that the detected sensitivity of transport characteristics of the considered nanoparticle network to 
TNT may be used to design a sensing nanodevice. 
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Molecular-scale conductors attract significant interest 
and attention in both fundamental and application re- 
search during the last two decades [H, [3, y, ll] ■ The un- 
ceasing interest of the research community to these con- 
ductors originates from their actual and potential useful- 
ness in inventing and designing various nanodevices. fn 
particular, the integration of nanotechnology with molec- 
ular recognition poses numerouspossibilities for creation 
of novel sensing devices 0, H, Significant 

efforts are being applied to study electron transport in 
networks made out of narioparticles connected by linker 

, 13, [3] ■ It is both interesting 



molecules 11, 12, 13, 14 




and important to investigate their potential in designing 
molecule recognizing sensors. 

The present work is motivated by recent experi- 
ments carried out to explore transport properties of a 
three-dimensional molecular network consisting of gold 
nanoparticles linked by oligophenylenevinulene (OPV) 
molecules. The molecular network was built using a cow- 
pea mosaic virus (CPMV) as a scaffold. The CPMV is 
a particle of an icosahedral geometry made out of sixty 
copies of a protein subunit whose average diameter is 
about 60nm, as was determined by X -ray crystallo- 
graphical analysis 17[. To use the CPMV as a scaffold 



for the molecular network it was genetically engineered 
in order to provide cystein residues at selected positions. 
Gold nanoparticles with diameters ~ 5nm were bound 
to these residues. Due to the icosahedral symmetry of 
the virus the residues and attached gold particles are ar- 
ranged into complex three dimensional structure of the 
same symmetry. The latter determines interparticle dis- 
tances. The CPMVs decorated with gold nanoparticles 
were exposed to OPV molecules which linked the gold 
particles producing conductive network on the viruses. 

To measure the conductance of a network built on a 
single virus the latter was assembled between proximal 
probe electrodes. The electrodes were made out of pair 
of lOOnrn wide gold leads patterned on a silicon wafer. 
After measuring the baseline conduction, the entire chip 
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FIG. 1: Color online: The nanoparticle network element and 
current-voltage characteristics for the network alone (upper 
blue line) and the network exposed to TNT vapors (lower red 
line) obtained in the experiments. 



was exposed to the target, and the transport measure- 
ments were repeated after binding the target. Actually 
the test chip included several pairs of electrodes which 
allowed statistical analysis of the results. Strong changes 
in the network conduction were observed when the lat- 
ter was exposed to trinitrotoluene (TNT) vapors. The 
exposure resulted in a significant decrease in the con- 
duction (Fig. 1). This conduction decrease occurred at 
rather small values of the bias voltage applied across the 
network (— 0.2y-^0.2y). The experimentally disclosed 
sensitivity of the above described molecular network to 
TNT makes it potentially useful to design a TNT sensing 
nanodevice. 

To theoretically analyze and explain these experimen- 
tal results we computed eigenenergies and correspond- 
ing wave functions for the molecular network element in- 
cluding two gold nanoparticles connected by the OPV 
molecule both with and without attached TNT. We 
did perform full self-consistent computations using the 
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NRLMOL software package [18[. In addition to self- 
consistently solving the Kohn-Sham equations we have 
optimized the geometries of the examined molecules. 
Once we arrived at the resulting Hamiltonian matrices we 
proceeded using Lowden's method [19] of symmetrical or- 
thogonormalization to build up atom centered Wannier- 
like wave functions .20] starting from the nonorthogonal 
gaussian orbitals. Using these as a new basis functions we 
reconstructed the Hamiltonian matrices. Now, the ma- 
trices became nearly diagonal enabling us to separate out 
three diagonalized blocks corresponding to the left and 
right gold nanoparticles {L, R) and the molecule (M) 
in between. Small off-diagonal elements in the resulting 
Hamiltonian matrices between the L — M and M — R 
blocks describe the coupling of the molecular linker to 
the gold nanoclusters in the considered network. 

The final matrices used in the following transport cal- 
culations are obtained separating out those parts of the 
Hamiltonian matrices which could be placed in the win- 
dow around the equilibrium Fermi energies Ep. The 
window width is determined by the value of the bias 
voltage. So, the full Hamiltonian matrices are replaced 
by these reduced matrices of relatively small dimension- 
alities. This procedure is justified for only those states 
whose eigenenergies are rather close to the Fermi energy 
could actually contribute to the electron transport. 

In our transport calculations we employ the effective 
Hamiltonian for the molecule linking two adjacent metal 
nanoparticles. The latter is written in the usual form 



(see e.g. Refs. [211,122!): 
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HL+Hn 



(1) 



Here, the term Hm corresponds to the molecule, and 
E[l.r describe the coupling of the molecule to the gold 
nanoclusters. Omitting for a while the last two terms 
and keeping in mind that the orbitals included into the 
basis set are orthogonormalized, we may introduce the 
retarded Green's function for the single molecule Gq {E) 
0. The latter is defined by the matrix equation: 



[iE + i7^)i-HM]G(liE)^i 



(2) 



where / is the identity matrix. The parameter rj de- 
scribes the broadening of the molecule energy levels due 
to random motions of the surroundings at finite tempera- 
ture T and takes on values of the order of thermal energy 
kT. 

To proceed we employ the Dyson equation relating the 
retarded Green's function of the molecule coupled to the 
metal nanoclusters G^{E) to the Green's function of a 
single molecule. The equation reads [l^: 



G^iE)^[A{E)-±{E)]- 



(3) 



Within the chosen basis set of the wavefunctions the ma- 
trix A{E) — [Go'(-E)] is a diagonal matrix 



A,j{E) = {E + ii] - E,)6,^ 



(4) 



where Ei are the eigenenergies of the molecule. 

In the Eq. ([3]), the self-energy term character- 
izes the coupling of the molecule to the left and right 
nanoclusters: 



j:{e)^j:l{e) + ±r{e). 



(5) 



The matrix elements of self-energy corrections have the 
form Q: 



^ E 

k 



(6) 



Here, /3 e L,R, Wik^/s are, respectively, the couphng 
strengths between i -th molecule state and "fc" -th state 
on the left/right metallic cluster, ek,i3 are energy levels of 
the corresponding nanoclusters, and the parameters (Tk^fs 
describe the thermal broadening of the electron levels at 
the clusters. The summation over k in the Eq. ^ is 
carried out over the states of the left/right cluster. 

When the bias voltage V is applied across the system 
shown in the Fig. 1, this causes charge redistribution, 
and subsequent changes in the energies Ei and ek,0- In 
consequence, the matrix elements Aij and values 
vary as V changes. This affects the electron transmis- 
sion function T given by the expression: 



(7) 



where f-^'^ = — 2ImEi_/f, and G^ is the advanced 
Green's function of the molecule (G"^ = (G^)^) . Never- 
theless, our computations showed that at low bias voltage 
values {\V\ < 0.251^) its effect on the electron transmis- 
sion is weak and may be neglected. Therefore, in further 
analysis we use the expression for the electron transmis- 
sion function computed at = 0. Wc remark, that the 
self-energy parts given by Eq. ([6]) depend on the values 
of Ck^fS representing electron structures of the gold nan- 
oclusters, and on temperature inserted in the terms ak,i3 ■ 
So, the changes in the electron structure of the clusters 
may affect the transmission as well as the changes in the 
electron structure of the linking molecule. The latter 
is included in the expressions for the Green's functions 
G«'^. 

Now, we employ the standard expression for the elec- 
tron tunneling current flowing through the molecule Q: 



7=4/ dET{E) [f{E ~ ^iL) - f{E - m)] . 



(8) 



Here, e is the electron charge, Ti is the Planck's con- 
stant, T[E) is the electron transmission given by Eq. 
([7|), f{E) is the Fermi distribution function for the en- 
ergy E. Chemical potentials ^ll^r are attributed to the 
left /right nanoparticles, and they are shifted from the 
equilibrium Fermi energy Ep due to the bias voltage V 
applied across the system: 



/ii = + (1 - v)eV, fiR = Ef 



veV 



(9) 
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FIG. 2: Color online: Current-voltage characteristics calcu- 
lated for a chain of nine gold nanoparticles linked by OPVs 
alone (solid lines) and by OPV-TNT complexes (dashed 
lines). The number of linkers N = 8. 



where v is the division parameter which shows how the 
voltage V is distributed between the nanoclusters. 

We apphed Eqs. (El)-® compute to the tunneUng 
electron current through a part of the experimentally 
studied network. The chosen part included eight pieces 
connected in series, each piece consisting of two gold 
nanoparticles linked by a molecule. The number of pieces 
roughly corresponds to the number of gold nanoparticles 
set along the half of the CPMV circumference. We as- 
sumed the bias voltage to be applied across the whole 
chain, so the voltage across two adjacent gold clusters 
appeared to be smaller than the net voltage V. Sup- 
posing that all molecular links in the chain are the OPV 
molecules and none of these are loaded with TNT, we did 
obtain the current-voltage characteristics shown as solid 
lines in the Fig. 2. At low bias voltage the characteristic 
is an ohmic curve which strongly resembles the experi- 
mental I — V curve (see Fig. 1). Both experimental and 
calculated I — V curves are symmetrical, and the com- 
puted current takes on values close to those observed at 
the same bias voltage, as follows from comparison of the 
figures 1 and 2 (left panel). When TNT molecules are 
attached to all OPVs included in the chain, our calcula- 
tions result in the critical reduction of the conductance at 
low bias voltage as presented in the Fig. 2. This agrees 
with the experimental results. 

We remark, that computed drop in conductance sig- 
nificantly exceeds that observed in the experiments. The 
mismatch could be reduced if we assume that TNT 
molecules are attached to some OPV linkers in the net- 
work not to everyone. Within the chosen model we can 
show that the decrease in the number of TNTs attached 
to the OPV linkers brings enhancement in the conduc- 
tance of the chain which manifests itself as the increase 
in the slope of the relevant I — V characteristics at low 
bias. The corresponding lines would be positioned be- 



tween the dashed and the solid lines in the left panel of 
the figure 2, and the smaller is the number of TNT-OPV 
complexes the closer is the I — V characteristics to the 
solid line. The latter corresponds to the chain of gold 
nanoparticles linked by OPVs alone. 

At higher voltage the conductance of the OPV-TNT 
molecular complexes enhances and becomes greater than 
the conductance of single OPV molecules. As fol- 
lows from our calculations, at the energies close to Ep 
self-energy terms S^.i? accept smaller values for joint 
molecular linkers than for single OPVs connecting gold 
nanoparticles in the network. The current between the 
nanoparticlies is defined by Eq. ([8|), and at low values 
of the applied voltage energies close to Ep give the pre- 
dominating contribution to the relevant integral. So, the 
reduction in S/,^ results in the reduction of the conduc- 
tance at low bias voltage. 

Finally, in the present work we studied transport prop- 
erties of the network consisting of gold nanoparticles 
linked by OPV molecules. These studies were motivated 
by the experimental results which revealed strong sensi- 
tivity of this system to TNT vapors. We carried out ab 
initio electronic structure calculations for the network 
element including two gold nanoclusters connected by 
the OPV molecule alone and by OPV-TNT molecules. 
Then we applied the results to calculate electron cur- 
rent through a chain made out of gold nanoclusters con- 
nected by molecular links. We remark that the simpli- 
fied structure of the network adopted here, prevented 
us form quantitative comparison of the computed re- 
sults with those experimentally obtained for the realistic 
three-dimensional network of rather complicated struc- 
ture. Nevertheless, we showed that at low bias volt- 
age the conductance of the considered chain strongly de- 
creases when TNT molecule becomes coupled to OPV 
linkers. This agrees with the experimental results. In 
whole, the present results reasonably explain the ob- 
served sensitivity of the above described nanoparticle 
networks to TNT. Also, they elucidate the potential of 
such networks in designing chemoselective sensing nan- 
odevices. 
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